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We present a hierarchy of discrete systems whose first members are the lattice modified 
■ Korteweg-de Vries equation, and the lattice modified Boussinesq equation. The iV-th member in 

the hierarchy is an iV-component system defined on an elementary plaquette in the 2-dimensional 
lattice. The system is multidimensionally consistent and a Lagrangian which respects this feature, 
i.e., which has the desirable closure property, is obtained. 
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Abstract 



^ ■ 1 Introduction 

What we now call the lattice modified Korteweg-de Vries (MKdV) equation appeared, up to a 
point-transformation, as early as 1894 in the work of Bianchi [4] as the permutability condition for 
Backlund transformations of the sine-Gordon equation, ft appeared again as an integrable partial 
difference equation via Hirota's method in [10] , and it was derived via the direct linearization 
scheme in |17| . 

To describe the lattice MKdV in the modern setting, let v = v(l, m) be the dependent variable 
depending on position (I, ra) in a 2-dimensional lattice, where I and m shift by integer values, 
jy-^ ' and let p, q be parameters associated with the I, m-lattice directions respectively. Then the lattice 

\ MKdV equation reads 

• . p(vv — vv) = q(vv — vv), (1) 

where " " denotes a shift in the Z-direction, and """ denotes a shift in the m-direction. 
£Nj ' A generalized direct linearization scheme proposed in |19| led to a higher order analogue of the 

lattice MKdV equation, which can be identified as a lattice version of the modified Boussinesq 
(MBSQ) equation of [21]. This is a scalar equation on a stencil of 9 lattice points, given by 



2 ~ 2 ~ \ ~ / 2 ~ 2 

p v — q v \ v / p v — q v \ v I v v \ I v v 



j^' V pv — qv J v \ pv — qv J v ~ \v | J ~ \v ^ 



= f "I "« " (2) 



A more convenient form of this equation appeared in [20] . specifically a coupled system of equations 
in two dependent variables v,w defined on an elementary plaquette: 

pv — qv pvw — qvw pw — qw . 

Z = ~ = ■ 

v vw w 

The scalar equation |2]) can be derived from ([3]) by eliminating one or other of the variables. 
As shown in [20], in either form this is a discrete analogue of the potential modified Boussinesq 
equation [2T] . Symmetries and conservation laws for the two-component form of the equation were 
derived in |22] , 

The lattice Gel'fand-Dikii (GD) hierarchy.that is the system given in [19] which goes to the 
continuous GD hierarchy under suitable limito is a multicomponent system whose first members 



the 1970s members of the Russian school, including Gel'fand and Dikii, studied higher order spectral problems 
generalizing aspects of the Korteweg-de Vries equation (c.f. e.g. [71ll6|); hence these types of systems have been referred 
to as Gel'fand-Dikii hierarchies. 
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are the lattice potential KdV equation and a lattice version of the Boussinesq equation. A lattice 
'modified' GD hierarchy was also proposed in the Appendix of [19], the first two members of which, 
((J} and p|. were given explicitly, but its higher members were left implicit due to the difficulty of 
writing corresponding scalar equations in explicit form on increasing stencil size. 

In this paper we study a hierarchy of systems whose first two members are (JXJ) and (J3j> - Different 
to the direct linearization approach of [19] it will be obtained by a reduction procedure from the 
discrete Kadomtsev-Petviashvili (KP) equation of Hirota [9], and in fact the precise relationship 
between the system studied here and those of [19] is not yet known beyond the first two cases. Like 
the first cases ([T]) and ©, each system is defined on an elementary plaquette of the 2-dimensional 
lattice (i.e., it is a quad equation); in particular each is consistent in multidimensions and possesses 
a Lax representation. Our main objectives are to present the derivation of this hierarchy as well 
as a Lagrangian formulation that respects its consistency property. 

The new kind of variational principle for multidimensionally consistent systems proposed in 
[13j relies on a closure property of the associated Lagrangian, allowing it to be interpreted as a 
closed form. A Lagrangian with this property was given in [15] for the generic member of the 
lattice GD hierarchy. We give here a Lagrangian for the generic member of the hierarchy in this 
paper, for which there is a clear connection to that of the discrete bilinear KP equation, whose 
Lagrangian was given in |14j . 

In Section [2] we derive the iV-component lattice system, give the Lax pair, describe its mul- 
tidimensional consistency and give some important symmetries. In Section [3] we present the La- 
grangian structure, and in Section [3] we demonstrate that this Lagrangian obeys a closure relation 
on solutions of the system. Section [5] contains some additional remarks. 



2 The system 

It was observed in [2] that the one-variable and two- variable equations ([l| and (J3| can be obtained 
by a natural reduction procedure from the discrete KP equation of Hirota [9]. Generalization of 
this reduction procedure then led to a quadrilateral lattice equation in N variables vi, . . . ,vn that 
can be written as 

pV n -!V n - qVn-XVn pv n v n+ i - qv„v n +i _ r-, , n ,.n 

5 = 5 > n £ {1, . . . ,N}, (4) 

V„-1V„ V n V n +l 

where by conventiorjf] Vq = fjv+i = 1- This system, which contains ([T]) and © as cases N — 1 and 
N = 2, is the main object studied in the present paper. 



2.1 Reduction from Hirota's discrete KP 

We begin by giving the reduction procedure from which Q may be obtained. The starting point 
is the Hirota discrete KP equation [9] 

f t = pff — qff, (5) 

where r = r(Z,m,n), f = t(1 + 1, m, n), f = t(1, m + 1, n), f = t(1, m, n + 1) etc. are values 
of the dependent variable r as a function of independent variables l,m,n £ Z. In this equation 
p,q 6 C\{0} are parameters, they are often taken equal to 1, which up to a non-autonomous change 
of variables is without loss of generality. Such change of variables will however not be possible after 
the reduction procedure. From the consistency property of ^ there can be obtained a Backlund 
transformation as follows: 

fa = pro — Xfa, (6a) 
r<7 = qfa — Xfa. (6b) 

2 This convention is a convenient notational device which avoids the need to separate equations such as {4} into sub 
cases n = 1, n £ {2, . . . , N — 1} and n = N. 
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Here, if r = r(i,m,n) is a solution of (|5j then the coupled system ([6]) for a is consistent and 
determines a function a — a(l, m, n) also satisfying (0). This is a weak Backlund transformation^, 
in fact the equation emerging as compatibility constraint for ([6]) is 

pff — qff 
ff 

and it can be directly verified that ((BJ is a Backlund transformation, in the strong sense, for 
this equation. The equation ([7]) lies on a 10-point stencil on the three-dimensional lattice, by 
integration it can be seen as a non-autonomous version of 

The equation ((4| governs solutions of (J7J satisfying the following additional constraints 

t(1, m, n + N + 1) = r(i, m, n), l,m,n£Z, (8a) 

r(Z,m,0) = l, l,m€Z, (8b) 

which can be seen by identifying 

Vn(l,m) := r(/,m,n), l,m£Z, n 6 {1, . . . , N}. (9) 

The first constraint (|8a[) is simply an imposed periodicity, the second constraint (|8b[l is not so 
usual, it consists of some particular boundary data. 



prr — qrr 



(7) 



2.2 Integrability 

The original Backlund transformation ((6} of equation (J7J) does not preserve the constraints (0 but 
it is consistent with the first constraint (|8afl . That is, we may write 

a(l, m, n + N + 1) = a(l, m, n), l,m,n£Z (10) 

without breaking the consistency. The further observation that this Backlund transformation is 
already linear in a makes it, through this reduction procedure, a Lax pair for Q. In other words 
defining 

(j>n(l,m) := a(l,m,n), !,m 6 Z, n 6 {0, . . . , N}, 
and (|10p reduces (O to the system 



and imposing 

where 
and 



<E> = L p (v, ■&)<!>, $ = Lq(u,u)$, 



$ := 



/ pvi/vx -\/vi 



L p (v,v) = 







(11) 
(12) 
(13) 



PV2/V2 —\vi/v 2 







V 



-Xv N - 1 /vn 
P 



(14) 



• • ■ • ■ ■ pvn/vn 

-\v N •■• 

is the (TV + l)-square Lax matrix. The consistency condition of (|12p 

L p (v, v)L q (v, v) — L q (v, v)L p (v, v) = (15) 

gives the system Q. 

The multidimensional consistency of (|4]), like the Lax pair, is basically inherited from ([5}- For 
this property, the vector v = v(l, m) on 1? is considered instead as a function on the d-dimensional 
regular lattice v — v(h, h, ■ ■ ■ , Id) satisfying system @ in all pairs of directions, where parameters 
p, q appearing in ((4]) associated with the original I and m directions are to be replaced with pairs 
of parameters taken from the larger set pi,...,pa associated with the corresponding h,...,ld 
directions. The consistency of this system can be verified directly by calculation. 



3 This terminology was probably first used by McCarthy [B] to distinguish the situation where an equation is strictly 
sufficient for compatibility of the Backlund system, as opposed to a strong Backlund transformation, where an equation 
emerges as both necessary and sufficient. 



3 



2.3 Point symmetries 

The system Q is invariant under the scaling symmetry 

(yi,...,vtr) — >• (m Unvn), m, ■ ■ ■ , G C \ {0}, (16) 

as well as the less obvious transformation 

(«i, V2, ■ ■ ■ ,vn-i,vn) — > («a/wi, V3/V1, . . . , vn/vi,1/vi). (17) 

The latter transformation here generates a group of N similar symmetries of Q. Such autonomous 
point symmetries allow an equation to take on alternative guises via a non-autonomous point trans- 
formation that preserves the autonomous nature of the equation, and also preserves its multidi- 
mensional consistency property (see for instance [3]). Explicitly, if Si and S2 are (commuting) 
transformations of the vector v — (vi, . . . ,vn) taken from the group generated by (|16[) and (|17|) . 
and v is governed by Q, then the equation governing v' = [S\ ■ S , ™](«) is both autonomous and 
multidimensionally consistent. The autonomous nature of the equation is preserved due to the 
fact the Si are symmetries, whilst the natural extension to multidimensions of the transformation: 
[S{ ■ S™] — ¥ [S l i ■ S l 2 ■ S3 3 •...], demonstrates also that the consistency is preserved. 

This is an important part of the transformation theory of multi-component systems such as Q 
because invariants of such transformations are not presently known beyond the multi-afEne class 
[I]. The less obvious symmetry (|17[) is included here because to distinguish between genuinely 
different quad equations in the multi-component class, it is necessary to be aware of all such 
transformations. 



3 Lagrangian structure 

For discrete systems where each lattice direction is on an equal footing it appears that the La- 
grangian is a more natural object than the Hamiltonian, given that it does not single out any 
one direction. Therefore an important feature of the multicomponent quad equation @ is its 
Lagrangian structure, which as we will see is inherited from that of Hirota's discrete KP equation 
([5]). The more remarkable feature is its natural interplay with the multidimensional consistency 
of the model, which will be explained in the following section (Section 2]) . 
A Lagrangian for the system ([4]) is 



qViV i+ i\ 1 / / qviv i+1 \\ 1, _ „ v 

~^-z I + — I In I — — + - (tavi mvi+i - hxvi lnu i+ i) 

pviVi+i J 4 \ \ pviv i+1 J J 2 



which through the discrete Euler-Lagrange equations gives 2 copies of the equations. The function 
Li2 is the dilogarithm function 

Li a (*) = - f Hl ~ z) dz, (19) 



which appears in many areas of physics, such as quantum electrodynamics (e.g. vacuum polar- 
ization) and electrical network problems; it also comes up in connection with algebraic K-theory, 
representation theory of infinite dimensional algebras, and combinatorics |llj . 
By discrete Euler-Lagrange equations we mean in this case 

^i = ^i +T -i^i +r -i^i = i€{1i ... )iV }, (20) 

where T v is a shift operator in the lattice direction associated with parameter p, and similarly T q 
is a shift operator in the lattice direction associated with parameter q. Then Tp ,Tq are simply 
backwards shifts in the lattice. 
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The discrete Euler-Lagrange equations arise from extremizing the action 



* - EEE{".(^) + K to (-^))' + i< h * to *«- , "* b '*«> 

+ KT)KI)}' ™ 

where the first two sums over I and m are understood to be taken over the whole 2-dimensional 
lattice. In fact, rather than the Lagrangian by itself, it is more instructive to look at the action, 
since it is on this level that we can easily see we are dealing with a reduction of the bilinear discrete 
KP equation 

A pq TT -(- A qr TT -\- A rp TT — 0; (22) 

here the dependent variable is r, and the A pq are antisymmetric constants (so that A qp — — A pq ) 
associated with the p, q lattice directions. Making the identification A pr /A pq = p, A qr /A pq = q 
and 7" — it, so that the equation is no longer covariant due to the singled-out lattice direction, 
brings us to equation ©. A Lagrangian for (I22C was presented in |14| . which, adapted for the 
non-covariant equation (JSJ), is equivalent to 

= EEEH$) + K'"(-SJ)) 2 + > f "-'"^> 

+ ln(l)ln(l)}. (23) 

It is clear that identifying r = Vi, f = Vi+i, and restricting the sum over n, we get the action ()21[) 
of the reduction. 



4 Closure property 

The Lagrangians for many integrable systems can be interpreted as closed forms; they obey a 
closure relation on solutions when consistently embedded in a higher-dimensional space. In a 
sense one can consider this to be multidimensional consistency on the level of the Lagrangian. It 
allows for a variational principle which includes the geometry of the space of independent variables, 
as well as that of the dependent variables; seeking extrema of the action under variations with 
respect to all of the variables, dependent and independent, leads to constraints which specify not 
only the equations of motion, but also to an extent the Lagrangians themselves. 

In the case of 2-dimensional systems, we have Lagrangian 2-forms, for instance Lagrangian 
functions which are evaluated on 2-dimensional plaquettes. For multidimensionally consistent 
systems we can suppose the action to be defined on a surface in an arbitrary number of dimensions, 
not just two, and this action will be a sum of the Lagrangian function evaluated on every plaquette 
in that surface. Demanding zero variation with respect to the surface implies a closure relation on 
the Lagrangian, while variation with respect to the dependent variable gives equations of motion. 

Examples can be found in both discrete and continuous cases of 1-dimensional systems [241 125] . 
a large class of 2-dimensional systems |13l 1151 23j , and also the 3-dimensional discrete bilinear KP 
equation |14| . In [15] it was demonstrated that the Lagrangian for an arbitrary member of the 
lattice GD hierarchy obeys the closure relation 

A p £ qr + A q C rp + A r C pq = 0, (24) 

on solutions to the lattice KP equation (the lattice KP equation is a consequence of the equations 
which constitute the lattice GD hierarchy itself). Here the operator A p is a difference operator in 
the p-direction of the lattice; if we again employ T p to denote the shift operator in the p-direction, 
and write the identity operator as id, then A p — T p — id. 
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To demonstrate the Lagrangian (|18[l satisfies the closure relation (|24|l . we need to invoke a 
5-term identity for the dilogarithm function 



Li 2 (s) + Li 2 (i)-Li 2 (st) = Li 2 ( 7 4)+Li 2 



st 



1 - st 

1-t 



st 



(25) 



which holds for s, t, st 7^ 1, up to imaginary constant terms, and also a 2-term identity which flips 
the argument 



Li 2 ( S )+Li 2 (Ij =-i(ln(- S )) 2 -^, 



(26) 



for s / 0. Further identities and other information about the dilogarithm function can be found 
in [12]. 

Observe that 



£pq ~f~ £qr 

N 



! t . ( qviVi +L v 



-Li 2 



pViVi+l 



+ L12 -rr-s + Li 2 



Li 2 



qviVi+i 



Li; 
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\rviv i+ i 



+5 I In 



+ i (ln 



1 n n [ gjA+i 
4 \ \ 



4 \ \ qviVi+i 



+ 2 (ln 



-4 ,ln 



+ i(lnii lnui + i - hauj lnt> i+ i) + i(lniji lnw i+ i - Iniii ln« i+ i) 

+ i(ln«i ln« i+ i - ln£>; lnu i+ i) - g (hi«i ln« i+ i - lnSj ln« i+ i) 

1 1 
— -(\nvi \nv i+ i - \nvi ln« i+ i) — -(in Si lnw i+ i - ln«, ln« i+ i) 



+ l n |£i±l N |lnf Vi 



ln(^hnf Vi 

Vi 
Vi_ 
Vi 



-ln|^±i]ln 



(27) 



Using the 2-term dilogarithm identity J26J) on the terms in the dashed-line boxes, and the 
5-point identity (|25p on the terms in the solid-line boxes, this becomes: 



G 



JV 

pvi qviVi+i — rviV i+1 \ I Vi+i rvtVi+i — pviVi+i 



~J I V r«i - qviVi+i / V v i+ i pviv i+ i — qviVi+i 



r« i+ i piw+i - gwm+i J \ Vi pviv i+1 - qvm+i 

v i+ i qiiV i+1 - rvSi+i \ ( qv~ t rv l v i +i - pViV i+ i 



+ ln( ^%+i yfi</i+i — ' ^%^%+i \ j / i/</t 



v i+ i pViV i+1 - qViV i+1 ) V m pv r v i+1 - qViV i+1 
! I _*i qVtVt+\ - rvm+i \^ n f qvi+i rvm+i - pvm+i 



Vi pviVi+i — qviVi+i J \ rvi+i pviV i+ i - qviV i+ i 
4 V V pviVi+iJJ 4\ V qviVi+iJJ 4\ V rviVi+i 



4 \ \ pviVi+i ) ) 4 \ \ (jWitii+i J J 4 \ \ rwm+i 
+ i(ln«i lnf^ + i - ln^i InSi+i) + i(lnui ln« i+ i - m«; lnu i+ i) 

+ i(lnw» ln« i+ i - lndi lnui+i) - i(lnii» lnv i+ i — In Si lnv i+1 ) 

— -(ln«i lnwi+i — Invi lnvj+i) — -(in Si ln-u;+i — lnwi InCi+i) 



From the equations, we have that 



pviVi+i - qViV i+ i vi 



(29a) 
(29b) 

f; - ' ' ■ — (29c) 

rvi — pvi v i+1 

for i £ {0, ...,N}, and of course vo — vn+i = 1 (so for i — 1 these are trivial identities). 
Identifying also vn+2 = wi, the equations hold trivially for i = iV + 1. Substituting these in gives 



pvi 


~ qvi 


Vi+1 


qviVi+i 


- rvm+i 


Vl 


qvi 


— rvi 


Vi+1 


rvtVi+i 


- pViV i+1 


Vl 


rvi 


- pvi 


Vt+1 


= Vn+1 


= 1 (so 


for i = 
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pvi+i qViVi+i — rvm + i \ I v i+1 rv l+1 v l+2 - pv l+l v l+2 



U. _ pvi+i gviVj+i — ■rvjVi+i \ ^ ^ 



rvi+i pviVi+i — qViVi+i J \ v i+ i pv i+1 v i+2 - qv i+ iVi +2 



rvi+i pviVi+i — qViVi+i J \ v t pviv i+1 — qviVi+i 

! , Vi+l qviVi+i - rViV i+1 N j/ qvi+i rv i+1 v i+2 - pv i+ iv i+2 ' 



v i+1 pViV i+1 - qViV i+1 J \ rv i+1 pv i+1 v i+2 - qv i+1 v i+2 _ 
ll _Vi_ qviVi+i - rVjVj+1 \,{ qi>i+i rViV i+ i - pViV i+1 



vi pViVi+i — qviVi+i / \ rvi+i pviVi+i — qviVi+i 



<i\ \ pv l v i+1 JJ 4\ V qiiv i+1 JJ 4\ V rvii i+1 

+ 1 L (_wA±A V _ 1 A ( _™m ±1 \ V _UJ psifii+i 



PViVi + 1 / / 4 V \ + 1 / / 4 V V + l 

+ i(ln ln-Si+i — ln^ilnDi+i) + i(lnvilnvi+i — lnvilnvi+i) 
+ i(lnfii In v i+ i — In ^ In Vi+i) — i (in Vi In v i+ i — In Vi In Vi+i) 
— i(ln ^ ln-Oi+i — lnfiilnui+i) — i (lnvi In Vj+i — lnvilnvi+i) 

+K^HtM*)KtMif)Kt)}- (30) 

Clearly some of these terms will cancel, reducing to 



N 1 



Li L > 



Vi+i rVi + iVi +2 —pv i+ iv i+2 
v i+1 pv i+1 v i+2 - qvi+\Vi +2 



Li2 



Vi rvjVi+i -pViVj+i 
Vi pViVi+x - qViV i+li 



+2 In 



lnp 



hi 



rvi — pvi 

pVx - qvi 
Vi+lVi+l 
Vi+lVi+l 



hi 



Vi+lVi+l 



hi 



Vi+lVi+l 
ViVi 
Vi Vi 



hi 



lnq 



I! 1 I', 



hi 



Vi+lVi+l 
Vi+lVi+l 



hi 



ViVi 
ViVi 



+ lnr[ln(Si + iUi + i{)i + i{)i + i) - In (v^iViVi)] - 2 lnr [ln(v i+1 v i+1 ) - ]n(yiVi)] 

+ [(hiv i+ i) 2 - (lnSi) 2 ] — [in ^i+i lnvi+i -lnvihiVi] + [lnSi+i In v i+ i — lnvihiVi] 

+2 [in Vi+i In 5i+i - In in In Vi] - 2 [in v i+ \ In v i+1 - In Vi In Vi] 

+2 [in i>i+i In Vi+i — In in In Vi] — 2 [in &+i In Sj+i — In Vi In Vi] 

— [in iii + i In Vi+i — lnvilnvi] + [lnUi+i lnvi+i — lnViln-Oi] + [(lnvi+1) 2 — (lnv 4 ) 2 ] 

— [in Vi+i In Vi+i — In Vi In «j] — [in Vj+i In v,:+i — In Vi In Vi] | . (31) 

On summation, all of these terms will cancel out, leaving zero. 



5 Further remarks 

Multi-component quad equations generalizing the one derived in [19] , as well as the one given here, 
are of course of interest. The multidimensional consistency was used as integrability criterion to 
obtain generalizations of the lattice Boussinesq equation in [H], and identification of these systems 
with a generalized dispersion relation within the direct linearization framework [26] opens up 
a potentially fruitful avenue in this direction. As indicated in Section 12.31 the transformation 
theory of such systems is quite involved even before non-local Backlund-type transformations are 
considered. 
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There is also a hierarchy presented in [2] with the lattice Schwarzian KdV and lattice Schwarzian 
Boussinesq equations as the lowest two members; this system was obtained directly in relation to 
the modified hierarchy (j4|) by generalizing a Backlund transformation known in the case N = 1 
since [15] : 

(32a) 



p v„v n 

1 V n+1 V n -i 

q v n v n 



(32b) 



where n g {1, . . . , iV} and by our convention no = vn+i — 1 as usual. It can readily be seen that 
this system is compatible in the z„ provided the v n satisfy Q, and similarly it is compatible in 
the v n provided the z n satisfy the system 



Zi+1 — Zi+1 Zi — Zi 



i 6 {!,... ,N- 1}, (33a) 



N+l ^ _ S 

g TT gj ~ gj _ ZN ~ ZN fqqM 

Unlike for the modified hierarchy from which it is derived, the expected Lagrangian formulation 
of this multidimensionally consistent Schwarzian hierarchy is presently not known. 

The Lagrangian formulation of the multicomponent quad equation in this paper contributes to 
the mounting evidence that the closure relation and the associated least-action principle are key 
properties of discrete integrable systems (and indeed of the corresponding continuous systems as 
well, cf. [23]). An important motivating factor to study such a structure is its potential physical 
significance; in particular the Lagrangian multi-form structure constitutes a possible departure 
point for the quantization of such integrable models along the line of a path integral formulation. 
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